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We investigate the one-triplet dispersion in a modified 
Shastry-Sutherland Model for SrCu2(B03)2 by means of a 
series expansion about the limit of strong dimerization. Our 
perturbative method is based on a continuous unitary trans- 
formation that maps the original Hamiltonian to an effective, 
energy quanta conserving block diagonal Hamiltonian TicS- 
The dispersion splits into two branches which are nearly de- 
generated. We analyse the symmetries of the model and show 
that space group operations are necessary to explain the de- 
generacy of the dispersion at k = and at the border of 
the magnetic Brillouin zone. Moreover, we investigate the 
behaviour of the dispersion for small |k| and compare our 
results to INS data. 

PACS numbers: 75.40. Gb, 75.50.Ee, 75.10.Jm 



I. INTRODUCTION 

The study of quantum spin systems exhibiting a fi- 
nite spin gap has advanced significantly through the re- 
cent synthesis of novel magnetic materials. While quasi- 
one dimensional systems have been studied for a long 
time now, new interest arose from the discovery of new 
quasi-two dimensional materials like CaV409 [Q and 
(V02)P207 %■ Particularly interesting is SrCua (603)2 
1^ since it is an experimental realization of the Shastry- 
Sutherland model HJH]. In Ref. B we presented an ex- 
tended version of that model (Fig. P and deduced its 
T = phase-diagram in the model parameter space. 




FIG. 1. An extract of the lattice we suggest for SrCu2(B03)2- 
The coupling J\ is assumed to be antiferromagnetic. The starting 
point of our analysis is the limit of strong dimerization ( J2, J3 — » 0). 



The Hamiltonian is given by 

H = Ji 2_^ Si ■ Sj +J2 2_^ Si ■ Sk +J3 2_^ Si ■ Si, 

<i,j> <i,k> <i,l> 



Ho 



Hi 



H2 



(1) 



where the bonds corresponding to interactions Ji , J2 and 
J3 are shown in Fig. |l|. For J3 = the model reduces to 
the original Shastry-Sutherland model. 
It can be easily verified that the singlet-dimer state (sin- 
glets on all strong bonds Ji, henceforth called dimers) 
is an exact eigen state of our model: a single spin in- 
teracts via J2 and J3 with 5 = objects only. Thus 
the corresponding terms in Hamiltonian (nl) do not con- 
tribute. The remaining expression simply gives — (3/8) Ji 
per spin. This article focuses on the region, where the 
singlet-dimer state is the ground state, called the dimer 
phase. In this phase a single excitation (magnon) is in- 
troduced by breaking up one singlet and substituting one 
triplet instead (the triplet's z-component is irrelevant). 
By hopping from dimer to dimer this triplet acquires a 
dispersion, which we intend to calculate. It suffices to do 
the calculations on an effective square lattice Toff, where 
one site represents one dimer. Note that Tes has a sublat- 
tice structure A/B, where A (B) corresponds to vertical 
(horizontal) dimers. 

We derive a power series expansion about the limit of 
strong dimerization {J2,J3 —^ 0) for the one-magnon 
dispersion. By fitting the expansion to INS-Data ob- 
tained by Kageyama et al. we deduce model parame- 
ters which show good agreement with parameters previ- 
ously determined. ||] 

The next section gives a short introduction to the method 
used. Before we calculate and discuss the one-magnon 
dispersion (Sec. IV) we use some of the model's symme- 



tries to derive interesting and useful relations between 
various hopping amplitudes (Sec. HI). 



II. METHOD AND QUALITATIVE PICTURES 

The one-magnon dispersion is calculated perturba- 
tively about the limit of isolated dimers using the flow 
equation method introduced previously P] . Given a per- 
turbation problem that can be formulated in the stan- 
dard way 



H = Ho+ xHs 



(2) 
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this method in its present formulation rehes only on two 
further prerequisites: 

(A) The unperturbed Hamiltonian Tig must have an 
equidistant spectrum bounded from below. The dif- 
ference between two successive levels is called an en- 
ergy quantum. 

(B) There is a number N 3 N > such that the per- 
turbing Hamiltonian Hg can be written as Ho, — 
'^n=-N -^" where T„ increments (or decrements, if 
n < 0) the number of energy quanta by n. 



The flow equation method maps the perturbed Hamil- 
tonian 7i by a continuous unitary transformation to an 
effective Hamiltonian Tieff, which conserves the number 
of energy quanta, i.e. [HcS^'Ho] = 0. Thus the effective 
Hamiltonian is block diagonal and has the form 



Tie 



Ti-o 



Y.^' E 



C{rn)T{m) 



(3) 



k=l 



\m\=k,M(m)=0 



where m is a vector of dimension k of which the compo- 
nents are in {±7V,±(7V-1), .. .±1, 0}; M(m) = signifies 
that the sum of the components vanishes which reflects 
the conservation of the number of energy quanta. The op- 
erator product T{ni) is defined by Tm — Tm-^Tm^ ■ • • T'm^ , 
where k is the order of the process. The coefficients C{m) 
are generally valid fractions, which we computed up to 
order fc = 15 for iV = 1 and up to order A: = 10 for N = 2 
(cf. Ref. §). 

We want to emphasize that the effective Hamiltonian Tictf 
with known coefhcients C{rn) can be used straightfor- 
wardly in all perturbative problems that meet conditions 
(A) and (B). The interested reader can find the C{rn) and 
additional information on our home-pages. For checking 
purposes we tested the C{m) by applying the method to 
toy models. In the iV = 1 case, for instance, we consid- 
ered the one dimensional harmonic oscillator perturbed 
by itself and an additional linear potential. Thus, with 
huj = 1 



Wo = i(P'+X2) = ata+l 



(4) 



ns = liP' 



X')+X = a^a + ^ + i=at + ^a . (5) 



To 



Ti 



In this case Eq. (0) can be solved exactly to give 



E„ 



a 



1 

"+2 



1 x^ 



21 



(6) 



We expand this last equation about a; = and obtain 

11 ^ °° 

£;„=n+- + -(2n + l)x--^(-l)V. (7) 



4 = 2 



By inserting the Ti defined in Eq. (|f) in the effective 
Hamiltonian (H) and calculating {n\Ti.cs\n) with ri e N 



we retain Eq. (|^) exactly up to IS'*' order. 

We now show that Hamiltonian (m meets conditions 
(A) and (B) for iV = 1. To this end we rewrite Eq. (0) 



f-» 



xHs , with 



Hs — Hi 



-H, 



•h 



and y 



Ji 



(8) 



(9) 



In the limit of isolated dimers {x = 0, with y/x finite) 
H is bounded from below and has an equidistant energy 
spectrum since up to a trivial constant Hq simply counts 
the number of excited dimers, i.e. energy quanta. 
To decompose Hs we follow the same proceeding as in 
Ref. PI and state the result 



H. 



r_i + To + Ti 



with 



^-^(i-DE'i.M 



To 



E^"M ' 



(10) 

(11) 
(12) 



where v denotes the pairs of adjacent dimers. For some 
fixed V we encounter the pair depicted in Fig. 0, the 
state of which is determined by \xi,X2)-, where a;i,a:;2 S 
{s^t^,t^ ,t~^} are singlets or one of the triplets occupying 
the vertical (horizontal) dimer, respectively. The super- 
script n G {0,±l}ini" stands for the S^ component. For 
the pair jxi, 0:2) we give the action of the local operators 
Ti in Tab. |. 



Tq I 


lt±\i±^) 


> 


|t±St±^) 


lt±\i") 


> 


jf«,i±l) 


|t±\i=Fl) 


—^ 


|iO,iO>-|t±\iTi) 


\t'\t°) 


> 


\t\t-^) + \t-\t^) 


\t°,t^') 


— > 


li±\i") 


^1 1 


|t±\.> 


> 


T|i",i±^)±|t±Si") 


\t'\s) 


— > 


|i\i-i)-|f-\ii) 



TABLE I. The action of the local operators Tq and Ti as they 
appear in Eqs. ( |ll|Jl2[ ) on all relevant states of the dimer pair de- 
picted in Fig. ti. These operators conserve the total S^ component. 
Note that Ti can only create another triplet on the horizontal dimer 
if there exists already one on the vertical dimer. This has also been 
noticed in Ref. [10]. Possible T±2 operators cancel out due to the 
lattice inherent frustration. Matrixelemonts not listed are zero. 



The remaining matrix elements can be constructed by 
using T2 = T-n- Note that we need to fix the orientation 
for singlets, say spin up with positive sign always at the 
right (upper) site of the dimers. Hence 7i and TLi ac- 
quire a global minus for oppositely oriented dimer pairs 
(reflection of the dimer pair in Fig. g about the vertical 
dimer). 

Instead of Hamiltonian (|l|) we can now use the effec- 
tive Hamiltonian Ticff (Eq- (m)) with the Ti defined in 
Eqs. (jl^jlj). The effective Hamiltonian (y) simplifies the 
computations considerably. Since [Tioffj^o] — 0; ^cff' is 
block diagonal allowing triplet conserving processes only. 
Thus the effective Hamiltonian acts in a much smaller 
Hilbert space than the original problem, which is a great 
advantage in the numerical implementation. 
In addition to this simplification the explicit form of 
Ticff' provides a simple and comprehensive picture of the 
physics involved. Imagine we were to put a triplet on one 
of the dimers in the lattice. In the real substance this lo- 
cal excitation would polarize its environment due to the 
exchange couplings and must be viewed as a dressed quasi 
particle surrounded by a cloud of virtual excitations fluc- 
tuating in space and time. The effective Hamiltonian in- 
cludes these fluctuations as virtual processes T{m), each 
weighted by the factor C{rn). Each process ends with a 
state having the same number of triplets as the initial 
state. Other quantum numbers such as the total spin 
are also conserved. As the order k is increased more and 
longer processes are allowed for and the accuracy of the 
results will be enhanced. Inspecting the weight factors 
C{m) shows that longer processes have less influence. 




FIG. 2. A pair of adjacent dimers (intra dimer coupling con- 
stant set to unity) connected by the perturbing interactions x and 
y. The local operators % as defined in Tab. Q acquire a global 
minus, if we reflect the pair on the indicated axis and keep the 
notation |zi,X2) as defined in the text. This is due to the singlet 
antisymmetry under refiection. 

Let US follow one of the possible virtual processes 
and understand why the observed triplet dispersion of 
SrCu2 (603)2 is rather flat (cf. Fig. ||). Suppose we be- 
gin with one triplet in the lattice as depicted in the upper 
left corner of Fig. |[ 




process ! 
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X 


\ - 


X 


— •— 1 1 
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FIG. 3. Leading (virtual) process for one-triplet hopping cor- 
responding to x^T{rn) = x'^T-iT-iT-iTiTiTi. Dark dots are 
triplets, bars are dimers. There is no lower order process leading 
to one-triplet motion. 

By applying Ti once we can create another triplet only 
on one of the two horizontally adjacent dimers as is clear 
from Tab. I and Fig. 0. From there we might create 
another one and so on till we can close a circle (bottom 
left state in Fig. ra). We now start destructing triplets 
by r_i processes and end up with the shifted triplet. 
The amplitude for this hopping is oc a:^. It is the largest 
amplitude one can find (see also Ref. |1^). Therefore the 
triplets are rather localized leading to a flat dispersion. 

Note that the Ti (T-i) processes are proportional to 
(x — y) (see Eq. (|ll|)). So, the leading triplet motion is 
0{{x — y)^) in leading order. Even local processes (with- 
out hopping) include at least two Ti (T-i) processes. 
Hence all perturbative amplitudes are at least of order 
(x — yY . We will make use of this fact later on. 

Examining the two-triplet sector S we showed that 
correlated hopping processes occur in second order al- 
ready. The actual dispersion, however, sets in only in 
third order. This much lower order [x"^ instead of x^) 
explains the much stronger two-magnon dispersion H . 

To quantify the picture constructed let |r) — |ri,r2) 
denote the state of the system with one triplet at r e T^s 
and singlets on all other sites. The amplitude t°^,_^ for a 
triplet-hopping from site r to site r' is given by 



Mr) 



(r'lT^efflr) 



(13) 



where the upper index o(r) € {u, h\ allows to distinguish 
whether the hopping started on a vertically oriented {v) 
or a horizontally oriented dimer {K). Furthermore we 
choose to split the hopping amplitudes into an average 
part is and an alternating part dts (s = r' — r) 



t: 



o(r) 



U + e'^'dt, 



(14) 



with Q — (tt, tt). 

The right hand side of Eq. ( |l3| ) can be easily implemented 
on a computer. For details see Ref. M . We want to point 
out, however, that all computations are done symboli- 
cally, i.e. we obtain results as functions (polynomials) of 
all model parameters. 



III. SYMMETRIES 



or simply 



Before we calculate the one-triplet dispersion quanti- 
tatively in the next section it is worthwhile to look at 
the symmetries the model in Fig. [^ displays. The two- 
dimensional space group of the model can be identified 
to be p4mm with the underlying point group 4inm as can 
be verified in Fig. H. 



o(r) 




U 2-fold axis 
H 4-fold axis 
mirror line 
glide line 

FIG. 4. A possible unit cell of the Shastry-Sutherland model 
with the dimers arranged at the sides of the cell. All symmetries 
are depicted. The two-dimensional crystallographic space group is 
p4mm. 

We will not treat all symmetry aspects but concentrate 
on those that will be of use later on. We choose the coor- 
dinate system parallel to the dimers such that one dimer 
(horizontal or vertical) lies in the origin and introduce 
two diagonals u and v crossing the origin with slope -I 
and I , respectively. The distance between the centers of 
two adjacent dimers is set to unity. 

Several relations between different hopping amplitudes 
can be deduced. To this end we define six symmetry op- 
erations which map the lattice onto itself. Note that the 
fixed singlet orientation can lead to negative phase fac- 
tors. 

ra^iy. Reflection about x/y-axis 

/: Inversion about the origin 

iV\-,T-il = I -r^.-Ti). 

o^ju- Reflection about the diagonals v/u plus translation 
by (0,-1) 

o-t,/„|ri,r2) = I ±r2,±ri - I). 

R: Rotation of 7r/2 about the origin plus translation by 

(0,-1) 

i?|ri,r2) = (-I)'-i+'-^|-r2,ri-I). 

Since these operations leave 7i unchanged they all com- 
mute with Ticff- Applying m^ to Eq. (n3) we find 



t°w = (-l)'^+'^t°J"" 

Analogously using niy, I, a^, (Jv and R we find 



(16) 



^o(r) ^ (_^).,+.,^o(0 ^^ 


(17) 


_ fO(r) 


(18) 


_ ,o(r-(0,l)) 
-S2.-S1 


(19) 


_ .o(r-(0,l)) 

S2,Si 


(20) 


_ / T\si+S2-i-O{r-(0,l)) 


(21) 



respectively. In particular Eqs. (jl^) and ( |18D yield to- 
gether 



o(r) _ 



(-si,0) 
o(r) 
(si,0) 
o(r) 

{0,S2) 



= (-i)"Ci: 



lj-0{r) _ ,o(r) 



(sl,0) 



(sl,0) 



0, if si odd, analogously 
0, if S2 odd. 



(22) 
(23) 



describing the interesting fact, that hopping along the 
axis has non zero amplitude only if this hopping moves 
the triplet an even number of sites (Foff). 



IV. DISPERSION 

Since Tics conserves the number of triplets the one- 
triplet dispersion is particularly easy to calculate. Start- 
ing with one triplet this excitation can only be shifted, 
i.e. the triplet hops on the effective lattice Foff. Addi- 
tional care has to be taken to account for the sublattice 
structure of Foff'. From Eq. (O) we get 



Hcff|r)=Ev^^l'- + r')- 



We introduce Fourier transformed states 

J{k+aQ)r 



■'>-7lP 



(24) 



(25) 



with the number of dimers L, the new quantum number 
a e {0,1} reflecting the sublattice structure and k a 
vector in the magnetic Brillouin zone (MBZ). Calculating 
the action of T-Ccs\(t, k) on these states yields 

Hcff k, k> = ^ ^(ir' + e''3"-dir')|r + r')e*(''+-«)'- 



(T,k) 



^dir'e-*(''+*«>' ^ 5] k + r')e^(''+*«)('^+'^') , (26) 



Vl 



|ff,k) 



with a — 1 — a. We used the definitions ( [1^ ) and (14). 
Further, since Eq. (Jl^) holds for even and odd ri + r2 
one has is — i-s and dig = di-s- Thus we can simpfify 
the sums over r' in Eq. ( Pq) to give 



Hofflcr, k) 



io + 



r>0 



cos((k + crQ)r) 



+ 



dto + 2^ dU cos((k + aQ)r) 

r>0 



|a,k), (27) 



with r > if and only if (ri > or ri = but r2 > 0). 
In Appendix A we show that dto — and dt^ — for 
ri + r2 odd. Hence 6 does not depend on a 



& = 2 ^ dtr cos(kr) 



(28) 



r>0 
ri-\-r2 even 



and Ticff is symmetric in the new states. The remaining 
2x2 problem can be solved easily to give the dispersion 

c^i/2(k) = ^^^ ±^v/(«o-ai)2+462 . (29) 

wo(k) 



start in 10*^ order only. One may object that b c onta ins 
dti^i, which could be larger, but from relation (A5) in 
Appendix A follows that dti^i=0. The amplitudes dt±2.o 
and dtp, ±2 in b also start in 10*'' order only. Hence the 
almost degeneracy in the one triplet sector can be under- 
stood on the basis of the symmetries of the lattice. 

By substituting y — in a-'o(k) we obtain the decimal 
numbers computed by Zheng et al. [|ll|| . Our series ex- 
pansion for the dispersion (E9) converges nicely but it 
can be improved further by the use of D-log Fade ap- 
proximants p2[ . 

Let us first consider the energy gap A :— w(0,0) as 
function of x and y. We fix the ratio of x and y and use 
D-log Fade approximants to extrapolate the remaining 
expression. For all x/y ratios we tested, we find the ap- 
proximants to be very stable, i.e. most of the possible 
approximants at a fixed ratio coincide very well. (Some 
are defective, i.e they show spurious singularities.) Fig. 
shows D-log Fade extrapolations for y = — a;tan(7r/6) 
and y = — a;tan(7r/8). In Ref. Q we used this technique 
to determine the line in the (x, ?/)-plane where the gap 
A vanishes. The vanishing of A indicates definitively the 
break down of the dimer phase. But it may happen that 
another excitation becomes soft before the elementary 
triplet vanishes (cf . Ref. js) ) or that a first order transi- 
tion takes place [§,|l^. In the present work we choose to 
examine the elementary triplet only. 



Thus the one-triplet dispersion splits into two branches. 
We want to point out, however, that at k = and at the 
borders of the MBZ (i.e. \kx + ky\ = tt or \ky — k^l = tt) 
the two branches fall onto each other leading to a two- 
fold degenerate dispersion. An analogous degeneracy is 
noticed in the two-triplet sector |^. In Appendix B we 
demonstrate that the degeneracy is due to the glide line 
symmetries R and cr„/„ and show that (ao — ai) and b 
both vanish. Moreover, (ao -I- ai) is a sum over r with 
ri+r2 even and thus contains "even" hopping amplitudes 
only. At k = and at the border of the MBZ we therefore 
find that one-triplet hopping takes place on one species 
of the two sublattices A/B in Fcff only. In other words, 
the triplets live either on the horizontal or on the vertical 
dimers. 

We calculated the amplitudes ir and dt^ (and therefore 
the dispersion) as exact polynomials in x and y up to and 
including 15*'' order. On appearance of this article these 
polynomials will be published in electronic form on our 
home pages. 

Expanding the square root in Eq. ( P9| ) about the limit of 
vanishing x and y produces terms oc x^y^ with a + f3 > 
10. Hence the energy splitting starts in 10*'' order and is 
negligible for all reasonable values of x and y. The fact 
that the splitting starts four orders later than the disper- 
sion may be understood by observing that (ap — ai) is a 
sum over r with ri -I- r2 odd. From the discussion at the 



end of section |IH| it is clear that t±J q and ^q ±i vanish so 



Mr) 




FIG. 5. D-log Pade extrapolations of the gap along two lines 
in the {x, y)-plane. Different approximants for a fixed ratio coincide 
very well. The {x,y) values where the gap vanishes constitute a line 
in the (x, j/)-plane indicating the definitive collapse of the dimer 
phase (cf. Fig. 2 in Ref. [6]). 

For y ~ — xtan(7r/6) the gap vanishes at x = 0.558(1) 
{y = —0.323(1)). At this point we get the lowest curve 
in Fig. |6| where we choose to show the dispersion along a 
triangle in the MBZ. 



that the leading "odd" process is t±2 ±i '^'' '•±1 ±2' which 
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FIG. 6. One-triplet dispersion in the MBZ at various points of 
the (x,j;)-plane. The lowest curve shows a;(k) at the point where 
the gap vanishes. With increasing y the dispersion decreases con- 
siderably as is clear from {u){'k, 0) — tij(0, 0)) oc (x — y) . 

In Sec. II we saw that all hopping amplitudes are pro- 
portional to {x — •yf' . Since the dispersion is a sum over 
these amplitudes it is clear that (a;(7r,0) — u;(0, 0)) oc 
[x — y)^. Indeed, for increasing y we see that the disper- 
sion decreases (Fig. ||). At a; = y we find the dispersion 
to be absolutely flat. This is a signiture of the fact that 
at X = y the total spin on each Ji bond is a conserved 
quantity [^ . Thus there will be no triplet motion. 

We turn to the behaviour of the dispersion a-'(k) for 
small |k| on the line in the {x, y)-plane where the gap van- 
ishes. After fixing the ratio x/y in Eq. (ES) and applying 
the D-log Pade technique we end up with an expression 



a;(a;;k) — exp 



Q(a;;k) 

/(a;;k) 



dx 



(30) 



where A is a function of k and P and Q are polynomi- 
als in X (where the leading coefficient is unity) of order 
M and N, respectively (shorthand: [M,N]-approximant). 
At k = the smallest positive zero of Q, say xq, is 
the point where the exponent diverges logarithmically to 
— oo, hence a;(a;o;0) = 0. We set ky = and will show 
that for small kr 



u>{xq; kx) oi k" + higher orders 
Differentiating Inw yields for kx —^ Q 



lim kx 



dk^f{x;kx)dx 



(31) 



(32) 



We decompose P and Q in linear factors 



P{x;kx) == {x - pi{kx)){x ~ p2{kx)) ■ ■ -{x-pMikx)) (33) 
Q{x; kx) ^ {x~ qi{kx)){x - q2{kx)) ■ ■ ■ {x - qNiK)) , (34) 

such that qi{kx)\k^=o = xq. In this way the factor 
(x — qi{kx)) dominates the behaviour of a for small kx- 
Further, since the dispersion is invariant under the sub- 
stitution kx —^ —kx we have A,pi,qi ex fc^ for small kx- 
Thus we are led to write qi{kx) — xq + pk"^, with positive 
j3. With these preparations we rewrite Eq. (p2[) 



a — lim kx 

k-,^0 



2(3kx 



{x-Xi)- Pkx) 



2^2 



g{x;kx)dx (35) 
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lim 2 



J^ g{xo - I3kl;kx) 

(l+y)2 



dy , 



(36) 



where we substituted x — xo—(5k'^y and g(x] k) for A{kx)- 
P{x; kx) ■ {x — qi)/Q{x] kx). Taking the limit yields 



2g(xo;0) 



1 



jdy = 2g{xo;0) 



(l + y)2 

A straightforward calculation shows that 

Pixo;0) 



gixo;0) = A{0) 



dxQ{xo,0) ' 



(37) 



(38) 



which we can easily calculate. Fig. uj shows the exponent 
a as function of the angle measured from the positive 
X-axis to the negative y-axis. 




(|)/3I 



0.50 



FIG. 7. The exponent a in u){kx) oc k" as function of the angle 
<j> between the x and (— y) axis. 

There is only one ratio x/y for which a — 1. At this 
point the dispersion vanishes with finite spin wave veloc- 
ity ((/> « 0.081(1) or a; « 0.679(1), y « -0.055(1)). 

Let us turn to a comparison of the theoretical disper- 
sion to experimental data for SrCu2(B03)2. To fit the 
dispersion to experimental data we make use of the pa- 
rameter dependance of our results by requiring the curves 



to go through certain points and solve the resulting set of 
equations. At k = (0,0) ESR 0, FIR (l| and INS § 
data suggest a value of ^(0,0) = 2.98nieV. At finite k 
we have to rely on the INS measurement, which contain 
rather large errors. In Fig. ra we show the INS data (bul- 
lets and error bars) and two of our fitted curves. For the 
parameter values given we get an excellent agreement. 



3.2 



2.9 



■ x=0.603, y=0, 



x=0.59, y=-0.()5. J,=6.67meV 



J,=6.16meV 




(0,0) 



(Jt,0) (it/2,7t/2) 



(0,0) 



FIG. 8. One-triplet dispersion. Our theoretical results fitted 
to INS data (bullets, errors at least as large as error bars). Due 
to the large errors it is not possible to fit the model parameters 
unambiguously. 

Because of the flatness of the dispersion and the com- 
parably large error bars it is not possible to fix the model 
parameters unambiguously . As sketched in Fig. g one 
can lower x and enlarge Ji and \y\ without loosing rea- 
sonable agreement. On the basis of the one-triplet disper- 
sion it is not possible to fix the model parameters quan- 
titatively. Our investigation in the two-triplet sector Iq] 
(where y is disregarded, since it is probably not very 
important) enables us to fix the parameters to the inter- 
vals {y ^ 0,) X = 0.603(3) and Ji = 6.16(10)meV. The 
corresponding one-triplet dispersion is the solid curve in 
Fig. |. 



V. CONCLUSION 

In summary, we have demonstrated the utility of our 
perturbation method introduced earlier |9[|. By the 
means of a continuous unitary transformation the orig- 
inal Hamiltonian is mapped onto a block diagonal, en- 
ergy quanta conserving, effective Hamiltonian Tieff- This 
allows to do the calculations in different energy sectors 
separately. 

Here we used this method to develop a picture of exci- 
tation processes in a model for SrCu2(B03)2 and calcu- 
lated the series expansion of the one-magnon dispersion 
uj about the limit of isolated dimers as a 15*'' order poly- 
nomial in the model parameters. 



We showed that the dispersion decomposes into two 
nearly degenerate branches with a splitting proportional 
to x^°. At k = and at the border of the magnetic 
Brillouin zone the two branches fall onto each other. By 
making use of our detailed analysis of the models symme- 
tries we showed that point group operations alone cannot 
explain these degenaracies. In fact we showed that the 
model's space group symmetries have to be taken into 
account. 

Moreover we analysed the critical behaviour of uj for 
small |k| at various ratios of x and y. It is found 
that the dispersion vanishes only for x ~ 0.679(1) and 
y « —0.055(1) with a finite spin velocity. 
Finally, we fitted our dispersion to INS data obtained by 
Kageyama et al. and found that it is not possible to 
fix the model parameters unambiguously from the one- 
triplet data alone due to very large experimental error 
bars. The parameter values determined by our investi- 
gations in the two-triplet sector ( g: y = 0, x = 0.603 , 
J = 6.67meV), however, lead to a one-magnon dispersion 
agreeing nicely with the experimental findings. 
The authors would like to thank Alexander Biihler for 
helpful discussions and the Regional Computing Centre 
of the University of Cologne (RRZK) and Dr. Kim Yong 
Talk in particular for generous assistance in our large 
scale computations. This work was supported by the 
DFG-Schwerpunkt 1073 and by the SFB 341. 



APPENDIX A: 

First we show that dtj.=0 if ri = r^. 
With Eq. (EQ) (the a^ symmetry) we have 



,o(r) _ ,o(r-(0,l)) 



Splitting both sides according to Eq. (14) we get 

leading to 

tg + dts = is2,si — dts2,si , for r even and 



Utc La 



dts2,si , for r odd. 



Taking the difference of both equations yields 



^Si,S2 — '^^S2,Si 5 



(Al) 



(A2) 



(A3) 
(A4) 



(A5) 



which proves the assertion. In particular dig = 0. 

We now show that dt^ = if ri -I- r2 is an odd number. 

According to Eq. (fiSh we have 



Mr) 



{v + r'\n,s\r) 
(r|Heff|r + r')-tt'+""'^ 



o(r+r') 



(A6) 



since Ti.cs has only real matrix elements in this basis (cf. 
Eq. d)). The last equality follows from Eq. ® . Split- 
ting both sides according to Eq. ( p!^ ) yields 



=> dtr' = e" 

=> dtr' = for r[ + r'^ odd. 



"''^"'dU, 



(A7) 



APPENDIX B: DEGENERACY OF w 

The dispersion (p9) will be two-fold degenerate if the 
square root vanishes. We will use the shorthand r 
even/odd for ri + r2 even/odd. 



k = 













A) b 
)rk = 


1 
= 

= Eq. ( 


2^ gives 








b 
2 


=y]du 

r>0 

r even 


ri=0 
r-2>0 

r even 


ri>0 
7-2=0 

r oven 


ri>0 
>-2>0 

r even 


ri>0 
r-2<0 

r even 



I1+I2 + I3 + h 



(Bl) 



We rewrite the first and the third sum on the right hand 
side of Eq. 



/l = 2^rfi-r2.ri and I^ — y ^dt^r2,rl 



(B2) 



ri>0 
T-2 = 



ri>0 
r9<0 

r even 



From Eq. (|2l|) (the i? symmetry) we deduce 



(J/Oip. ip^ — CtL — 'piy T'l 5 



and see that /i 
6 = 0. 



-I2 and /s 



-/4. Hence we find 



B) (ao - ai) = 
For k = we have (cf. Eq. (p^) 



ao - ai 



Y^ U [1 - cos(^(ri + ra))] = 2 ^ tr , (B3) 



r>0 



r>0 

odd 



and by recalling that dU = for r odd Eq. (|21| ) yields 



^ri,r2 ^— r2,ri 7 



(B4) 



so that we can use the same splitting as in Eq. ( |B1[ ) . Thus 
the energy-degeneracy at k = is due to the rotational 
symmetry R as defined in Sec III. 



II) kx + ky ^ n 

A) 6 = 
Making use of dtr = c?i-r, which follows from Eq. (18) 
we have for ky = n — k^ 



22 dtr cos{kx{ri - r2) + TTr2) 

V even 

yj dir2,ri cos(fc2,(r2 — ri)) cos(7rri) . (B5) 



In the last step we choose to rearrange the sum and ob- 
serve that if r is even we have ri and r2 both odd or both 
even. In both cases the identity 



cos(7rri) = cos(7rr2) 



(B6) 



holds. Inserting relation ( |A5D in the last row of Eq. (P35| ) 
we end up with 

- = - ^ (iirCOs(fc:r(ri -r2))cos(7rr2) , (B7) 



resulting in 6 = —6 and thus 6 = 0. 

B) (ao - ai) = 
Analogously to I) B) we have here 

= 2_, ^r [cos(kr) — cos(kr + 7r(ri + r2))] 

r 

= 2^ tr cos{kx{ri - r2) + 7rr2) 

r odd 

= 2^ ir2,ri cos{k^{r2 -ri))cos{7rri) 

r odd 

= -2^ iri,r2 cos(fc^(ri - r2)) cos(7rr2) , (B8) 

r odd 

where the last but one equality follows from Eq. ( pO| ) (the 
ay symmetry), i.e tri,r2 = tr2.ri, and from the fact that 
if r is odd we have that ri odd and r2 even or ri even 
and r2 odd. From that we see cos(7rri) = — cos(7rr2). 
The degeneracy over the remaining three borders of the 
magnetic Brillouin zone can be shown analogously. It 
is interesting to note, that the calculations necessarily 
involved glide line operations. The degeneracies can thus 
not be explained by considering point group symmetries 
only. 
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